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The in-plane resistivity, Hall resistivity and magnetization of graphite were investigated in 
pulsed magnetic fields applied along the c-axis. The Hall resistivity approaches zero at around 
53 T where the in-plane and out-of-plane resistivities steeply decrease. The differential mag¬ 
netization also shows an anomaly at around this field with a similar amplitude compared to 
that of de Haas-van Alphen oscillations at lower fields. This transition field appears insensi¬ 
tive to disorder, but reduces with doping holes. These results suggest the realization of the 
quantum limit states above 53 T. As a plausible explanation for the observed gapped out-of¬ 
plane conduction above 53 T, the emergence of the excitonic BCS-like state in graphite is 
proposed. 
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1. Introduction 

Graphite consists of stacked graphene layers along the c-axis, and is known as a com¬ 
pensated semimetal having almost the same number of electron and hole carriers. The carrier 
densities are about 3 x 10 18 cm~ 3 for both carriers, which is considerably smaller than that 
in typical metals. Elongated Fermi surfaces are located along the H-K-H (H’-K’-H’) points 
in the six edges of the hexagonal Brillouin zone. The electron-like Fermi surfaces exist near 
the K (K’) points and the hole-like ones near the H (H’) points. The band structure along 
the H-K-H (H’-K’-H’) points is well explained by the Slonczewski-Weiss-McClure (SWM) 
model. 1 ' 2 

In magnetic fields applied parallel to the c-axis, the band structure is quantized into Fan- 
dau sub-bands, which show dispersion only along the c-axis. Owing to the small effective 
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mass of carriers in the basal plane, the quasi-quantum limit state, where only the lowest 
electron-like (Landau index of n = 0) and hole-like in = -1) Landau sub-bands (each of 
them spin-split) cross the Fermi level, is realized in a moderate magnetic field of about 7.4 T. 

The application of higher magnetic fields to the quasi-quantum limit state in graphite 
causes an abrupt increase in the in-plane resistance above ~ 25 T. 3 Since the onset field 
strongly depends on temperature, this anomaly is ascribed to an electronic phase transition 
involving many-body effects. 4 A theoretical study by Yoshioka and Fukuyama (YF) sug¬ 
gested that the formation of a charge-density-wave state caused by Fermi-surface nesting in 
the n = 0 spin-up (Of) sub-band, and reasonably reproduced the temperature dependence 
of the transition field. 5 Here, the direct Coulomb interaction was eliminated by assuming 
the anti-phase charge modulation along the H-K-H and H’-K’-H’ edges. According to this 
model, the transition temperature for the density-wave state T c is given in the mean-field 
approximation as 


k#T c — 4.53 E i 


cos 2 (c 0 k F0T /2) 


exp 


N<n(E F )u(e)J’ 


( 1 ) 


COS(CokFOr) 

where E F is the Fermi energy measured from the bottom of the Landau sub-bands as shown in 
Fig. 1(a) and the No-f(E F ) is the density of states at the E v for the Of sub-band. k F o-f is the Fermi 
wavevector for this sub-band and c 0 is the lattice constant along the c-axis. u(e) represents the 
relevant pairing interaction as a function of the dielectric constant (e). In this formula, all the 
parameters E F , k F 0 j, M>t(£V), and u(e) depend on magnetic field. On the other hand, Takahashi 
and Takada claimed that the spin-density-wave state, characterized by another nesting vector, 
becomes stable when the spatial separation of electrons along the H-K-H and H’-K’-H edges 
is taken into account. 6 Since there is no report on the direct experimental evidence of the 
nesting vector, the real nature in this field region remains as an open question. 

While Eq. (1) is given for a particular 2k F -instability, we remark that Eq. (1) has a similar 
functional form to the general BCS-type formula for mean-field-type pairing 7 


fc B r c = L13E F exp(---), (2) 

N(E f )V) 

where N(E F ) is the density of states at the Fermi level and V is the relevant pairing interaction. 
When the Fermi energy approaches zero, T c goes to zero in Eq. (2). Consequently, a density- 
wave state with T c expressed by the Eq. (2) is expected to exhibit a reentrant transition back 
to the normal state if E F of the relevant sub-band approaches zero in higher magnetic fields. 
Takada and Goto calculated the renormalized band structure taking into account self-energy 
effects, and found that the 0T and -1 1 Landau subbands detouch the Fermi level almost 
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simultaneously at a crossing field of ~ 53 T s (see Fig. 1(b)). 

In fact, measurements of the in-plane resistance in magnetic fields of up to ~ 55 T ob¬ 
tained a clear indication of such a reentrant transition; the in-plane and out-of-plane resistance 
shows a sharp bend at a field close to the crossing field of ~ 53 T. 11,12 

Recently Fauque et al. investigated the longitudinal magneto-resistance (p zz ) in fields up 
to 80 T, and found an additional enhancement of the p zz above 53 T and successive reduc¬ 
tion of it above 75 T. 14 The authors interpreted these anomalies as sequential depopulation 
of Landau subbands, depopulation of hole-like -11 subband at 53 T and electron-like Of 
one at 75 T, and the emergence of novel density wave state in the intermediate field contrary 
to the preceding scenario. In this field region, the out-of-plane conduction exhibits gapped 
behavior, while the in-plane one (p xx ) remains metallic. This anomalous anisotropic conduc¬ 
tion was explained with assuming parallel conduction through the chiral edge state, which is 
characteristic of the quantum Hall state. 

To elucidate the true nature of the electronic states of graphite, studies of the Hall resis¬ 
tance in high magnetic fields will be important. In highly oriented pyrolytic graphite (HOPG), 
Kopelevich et al. observed multiple anomalies in Hall resistivity in field region between 16 
T and 43 T, and ascribed them to fractional quantum Hall plateaus. 15 They also reported 
the sign reversal in Hall resistance at 43 T, and interpreted as the change of majority carrier 
from electron to hole. Later, Kumar et al. reported that similar sign reversal occurs even at 
room temperature, 16 suggesting the irrelevance of the sign reversal of Hall resistance and the 
density wave state in HOPG. In experiments on single crystalline Kish graphite up to 30 T, 
the absolute value of Hall resistance steeply decreases in the density wave state. 17 The result 
does not show any feature of fractional Hall plateaus at least up to 30 T. The Hall resistance 
of single crystalline graphite above 30 T has not been reported as yet. 

To obtain further understanding of the true electronic states of graphite in the field region 
over 30 T, we carried out simultaneous measurements of the magneto-resistance and Hall 
resistance up to 65 T and of the magnetization up to 75 T. 

2. Experiment 

In this study, we investigated the transport and magnetic properties in various types of 
graphite in high magnetic fields applied along the c-axis. We investigated three batches of 
Kish graphite (KOI, K02, and K03), one super graphite (SG) made by a graphitization of 
piled polyimide sheets, and one natural New York graphite (NY). 

Since the in-plane chemical bonding is significantly stronger than the out-of-plane cou- 
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Fig. 1 . (Color online) Schematic drawings of the band structures (a) below and (b) above 53 T calculated by 
TG. S The horizontal solid lines indicate the Fermi level of pristine graphite while the dotted line in (b) indicates 
the Fermi level of the hole-doped one. SEM pictures of the sample edges (c) cut by razor and (d) as-grown crystal 
of KOI samples, (e) The out-of-plane magnetoresistance p zz of as-grown KOI, natural NY, and SG samples. 


pling, mechanically shaped graphite usually has disordered edges as shown in Fig. 1(c). 
Therefore, we mainly used as-grown crystals having clean edges as shown in Fig. 1(d) to 
extract intrinsic transport properties. In addition, graphite is known as a highly anisotropic 
conductor showing non-linear I-V characteristics in the field-induced density-wave state. 18,19 
To avoid possible mixing of in-plane and out-of-plane conductions, 20,21 we used thin samples 
with the ratio of length/thickness ~ 50 - 100 and formed electric contacts on the top surfaces 
of them for in-plane transport measurements. In this configuration, possible partial cleavage 
in the sample might affect in the evaluation of the absolute values of p xx and p xy . Thereby, we 
do not discuss quantitative aspect of the resistivity in this paper. 

Figure 1(e) shows the out-of-plane magneto-resistance ( p zz ) of KOI, NY, and SG. All the 
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samples show additional increase in p zz above 53 T as reported previously. 11,14 Therefore, we 
regard this enhancement as intrinsic nature in graphite. 

In-plane ( p xx ) and Hall resistance (p xy ) of the samples were measured simultaneously by a 
five-probe method using the numerical lock-in technique 22 at a frequency of 100 kHz or 200 
kHz. Epo-tek H20E silver epoxy and 30 //m-thick gold wires were used to form electrical 
contacts. Magnetization ( M ) was measured by an induction method using two pickup coils 
placed coaxially. To achieve large signal-to-noise ratio in the magnetization measurement, 
the crystals of KOI were piled up with their c-axis aligned as shown in the inset of Fig. 2(a). 
Magnetic fields were applied along the c-axis throughout the whole measurements. 

Pulsed high magnetic fields were generated using non-destructive-type magnets installed 
at The Institute for Solid State Physics at The University of Tokyo. For magnetotransport 
measurements, magnetic fields up to 56 T and 65 T (pulse duration of 36 ms and 4 ms, 
respectively) were generated using bipolar pulse magnets, which can generate both positive 
and negative fields. Mixing of the p xx component into the nominal p xy signal was eliminated 
by subtracting the data for different field polarities. Magnetization and some of p zz were 
measured up to 75 T using a unipolar pulse magnet with a duration of 4 ms. 

3. Results and discussion 

Figure 2(a) shows the field dependence of magnetization (M) of KOI. Monotonic sub- 
linear diamagnetism is observed up to 74 T. Figure 2(b) shows the differential magnetization 
dM/dH obtained by averaging seven data sets acquired under repeated experiments up to 62 
T. The dM/dH curve up to 11 T is shown in the inset of Fig. 2(b). Oscillations in dM/dH 
observed in the region up to 7.4 T corresponds to the de Haas-van Alphen (dHvA) effect 
indicating all sub-bands except for n = 0 and -1 detach from the Fermi level at fields below 
7.4 T. In the main panel, we can hardly identify anomaly at the transition field a shown below. 
On the other hand, a non-monotonic change was observed at around 53 T. 

Figure 3(a) shows the in-plane magneto-resistance p xx of KOI at various temperatures 
between 1.3 K and 4.2 K. Shubnikov-de Haas oscillations are observed in magnetic fields 
below 7.4 T. In higher fields, the observed abrupt increase and steep decrease in p xx (marked 
with a and a') have been regarded as the emergence and collapse of the density wave state, 
respectively. The critical field of the a transition shows strong temperature dependence in¬ 
dicating the importance of the many-body effect. On the other hand, the critical field of the 
a' transition show rather weak temperature dependence as reported previously. 11 Meanwhile, 
the other anomaly seen between a and a' pointed out in the previous reports 11,13 (which 
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Fig. 2. (Color online) (a) The field dependence of magnetization of KOI sample up to 74 T at 1.3 K. The 
inset shows a picture of the sample used for the magnetization measurements, (b) The field dependence of 
differential magnetization of KOI up to 62 T. a and a' indicate the fields where in-plane magneto-resistivity 
shows anomalies, while the open arrow indicates the field at which graphite goes into quasi quantum limit state. 
The inset shows the differential magnetization in the low field region. 


is referred to as /? anomaly below) is quite small in p xx . We observed clear ft anomaly for 
the relatively thick samples (not shown). This geometrical effect implies that the /3 anomaly 
might be the superposition of the out-of-plane resistivity. We will not go into details about 
the (3 anomaly in this study. The p xx above 53 T decreases with decreasing the temperature 
contrary to the gapped behavior in p zz , as reported previously. 14 

Figure 3(b) shows the field dependence of Hall resistivity p ^ at temperatures between 
1.3 K and 4.2 K. Shubnikov-de Haas oscillations are also observed in the below 7.4 T. 
In higher fields, an anomaly is clearly observed at a, but is less pronounced than that in p xx . 
Some readers might concern the relevance of this anomaly with the quantum Hall plateau. 
Since the normalized Hall resistance per single graphite layer at a varies from sample to 
sample in the other sets of experiments on KOI, we consider that this plateau-like anomaly 
at around p xx , -2 mf ! cm is not related to the quantum Hall plateau. Fractional plateau-like 
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Fig. 3. (Color online) (a) In-plane resistivity and (b) Hall resistivity of KOI sample at various temperatures 
as functions of magnetic fields applied along the c-axis. The inset of (b) shows the Hall resistivity up to 65 T at 
1.4 K. The result shown in inset was obtained for a different piece of KOI from that in the main panel. 


structures, which have been reported on HOPG, 15 were not observed in this measurement. We 
also observed that p^ gradually goes to zero as the applied field approached 53 T. As shown 
in the inset of Fig. 3(b), p xy remains almost zero even deeply in the high field state between 
53 T and 74 T. 


Figure 4 shows the field dependence of the conductivities cr xx and a xy of KOI. The cr xx and 
cr xy were obtained from the in-plane resistivity p xx and Hall resistivity p xy using the following 
relationship: 


CT vr = 


plx + ply 


CT 


Pxy 


xy 


pIx+pI ' 


( 3 ) 


( 4 ) 


'xy 


In high magnetic fields, the imbalance between the carrier densities of electron (n e ) and hole 
(n h ) is deduced within the framework of semi-classical simple two carrier model as follows: 


n h = 


0~xyB 


( 5 ) 


The field dependence of cr^B and thus estimated imbalance are shown in the inset of Fig. 
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Fig. 4. (Color online) The field dependence of conductivities cr xx (red) and cr xy (blue) of KOI sample at 1.3 
K. The inset shows the field dependence of cr xy B up to 30 T at 1.3 K. The scale on the right side of inset gives a 
difference between electron- and hole-density. 


4. The crxyB is nearly constant in the field region between 8 T and 15 T, then shows a slight 
decrease between 15 T and 25 T. This reduction has been ascribed to the effect of magnetic 
freeze-out, which diminishes carriers provided by ionized impurities. 24 As shown in the inset 
of Fig. 4, crxyB shows a dramatic decrease at a. As mentioned above, Uji el al. 11 ascribed this 
anomaly to the opening of a gap in the electron-like sub-band, which leads to the sign reversal 
of the Hall resistivity at higher fields. Contrary to this expectation, no apparent sign reversal 
was observed in our measurements up to 65 T. To clarify whether these features are intrinsic 
in graphite, we measured p xx and p xy in various types of graphite [Figs. 5(a), (b)]. In p xx , 
only the a transition field of SG is clearly higher than the other samples. The details of this 
difference will be discussed later. The p xy of KOI, K02, K03, and SG are negative, whereas 
only natural NY shows positive Hall resistance. Regardless of its sign, the p xy approaches 
zero at the a' transition in all the samples, suggesting it is intrinsic behavior in graphite. No 
quantum-Hall plateau-like structure was observed in any samples. 

Figure 5(c) shows the field dependence of the transition temperature for a ( T a ). According 
to the BCS-like formula in Eq. (2), and assuming the enhancement of the N(E F ) in proportion 
to magnetic field and ignoring field dependence of E v and V, Eq. (2) is rewritten in simpler 
form: 


T a = T* exp 



( 6 ) 


where T* and B* are constants. 4 This formula gives linear relation between In T a and the 
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Fig. 5. (Color online) (a) The in-plane magneto-resistance p xx and (b) Hall resistance p xy of five different 
samples at 1.3 K. (c) The field dependence of a transition temperature T a as a function of inversed field. The 
open circles and squares represent the reported results of Kish graphite irradiated with fast neutron for two 
(NI#2) and four (NI#4) hours, respectively. 9 (d) The field dependence of a' transition temperature T a > of five 
different pristine and two irradiated 10 graphite samples. NI#2 and NI#4 are irradiated with first neutron under 
the same condition, but different pieces from those in (c). The inset shows enlarged view of (a) at around the a' 
transition field. 


inversed field. Solid lines in Fig. 5(c) represent fitting curves for KOI, K02, K03, and NY by 
Eq. (6) with using parameters of T* and B* as listed in Table I. Here, the T* and B* simply 
modify the scales of vertical and horizontal axes, respectively. In our samples, only the SG 
exhibits non-linear phase boundary that cannot be pronounced by Eq. (6). Similar discrepancy 
was also observed in neutron-irradiated graphite. 9 Neutron irradiation to graphite is known 
to introduce lattice defects that act as acceptors, and hence, dope hole carriers. The reported 
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Table I. The list of samples investigated in this report. T*, B*. and r are fitting parameters in Eqs. (6) and (7). 


sample 

measured quantity 

r (K) 

B* (T) 

r (x 10- 12 s) 

paper 

KOI (Kish graphite) 

Pxx, Pxy, Pzz, M 

290 

144 

- 

this work 

K02 (Kish graphite) 

P XX 5 Pxy 

270 

136 

- 

this work 

K03 (Kish graphite) 

Pxxi Pxy 

257 

142 

- 

this work 

NY (natural graphite) 

Pxxi Pxy-9 Pzz 

192 

131 

- 

this work 

SG (super graphite) 

Pxx ^ Pxy > Pzz 

270 

136 

4.5 

this work 

NI#2 (Kish graphite) 

Pxx 

191 

128 

3.9 

Ref. 9 and 10 

NI#4 (Kish graphite) 

Pxx 

191 

128 

2.0 

Ref. 9 and 10 


field dependence of T a of neutron-irradiated Kish graphite are also shown in Fig. 5(c). 9 The 
results for irradiated samples significantly deviate from that of Eq. (6) similar to the cases 
of SG. Yaguchi et al. explained the qualitative difference of the transition temperature with 
introducing the pair breaking effect on density waves caused by charged impurities, 9 similar 
to the arguments on different types of graphite below 1 K. 23 In the mean-field form, the 
reduced transition temperature by the pair breaking effect (T 06 ) is expressed as, 23 


In 


y^PB 

1 a 




i 


n 


2 2 nk B T™T' 


(7) 


ijjix) is the digamma function, T„ is the transition temperature without the pair breaking effect, 
which is given by Eq. (6), and r is the relaxation time of pair breaking. The deviation from 
the functional form of Eq. (6) seen in neutron-irradiated Kish graphite is well reproduced by 
Eq. (7) as shown in Fig. 5(c). Apparently Eq. (7) also accounts for the reduction of T a seen 
in SG. With using the values of T* = 270 K and B* = 136 T taking from those in the K02, 
the phase boundary of SG can be reasonably reproduced with r = 4.5 x 10 -12 s as shown by 
the dotted line in Fig. 5(c). Although the microscopic origin is not clear, the T a in the SG is 
significantly affected by the pair breaking effect compared with the other four samples. 

Although the T a of SG shows the similar trace with those of neutron-irradiated samples, 
the transition temperature for a ’ (T a >) shows up significantly different. The T a > for all samples 
are plotted as a function of B in Fig. 5(d). The five pristine samples show almost the similar 
boundary, whereas the neutron-irradiated ones, quoted from literature, 10 show transitions at 
smaller fields. We consider that the reduction of the transition field at a' at a given temperature 
is not caused by the pair breaking effect because the observed reduction of that of SG is almost 
the same value as pristine Kish samples. The neutron-irradiated samples are known to have 
holes with the order of 10 18 cm -3 , which is hundred times lager than the change imbalance 
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for the pristine samples. Therefore, the reduction of the re-entrant transition field at a' in 
irradiated samples can be caused by the hole-doping rather than the effect of pair breaking. 

Let us discuss the transition at a' in the following. Firstly, since the differential magneti¬ 
zation shows anomaly at around 53 T in the same order of magnitude with those in the dHvA 
oscillations at low fields (open arrows in Fig. 2). We regard this anomaly as the feature of de¬ 
population of the Landau sub-band. Within the accuracy of the these measurements, however, 
we cannot determine the number of the depopulated sub-bands at ~ 53 T. 

To consider the number of populated Landau sub-bands above 53 T, we discuss the result 
of Hall resistance beyond the simple two-carrier model. For more precise argument, we have 
to calculate the in-plane and Hall conductivity using the Kubo formula: 


-i £ 


<x 




nkcr 
E'n < '^ F < '^«+1 


<77 1 j x 1 72 + 1) <72 + 1| j y \n) ~ <72 + 1 | j x | 72 ) <72 1 j y |77 + 1) 

(E n+] ~ E n f 


(8) 


In this calculation, one has to evaluate the summation of matrix elements for current operator 
j such as <72 1 j 1 72 + 1) or (n + 1| j \n), where |t2> or \n + 1) denote the states connected by this 
operator. In the present case of graphite in the quasi-quantum limit, these states should be 
chosen from the sub-bands having different Landau indices by number of 1, and same wave 
numbers and spins. Since the conductivity is determined by the summation of these matrix 
elements of the inter-sub-band transitions over a wide region of the reciprocal space, gap 
formation at a part of the sub-band may cause the limited contribution to the Hall conductivity 
in the density wave state while the gap is considerably smaller than the energy difference in 
the adjacent sub-bands. In other words, comparison between the experimental data of o- xy and 
such calculations will provide in-depth evaluation for several models of sub-band structures 
above 53 T. We leave it as an open problem because this calculation is beyond the scope of 
the present study. 

In the present study, we consider the possible sub-band structure through more qualitative 
argument. Firstly the observed small p xy seems to fit with the symmetric sub-band structure 
with respect to the Fermi level as illustrated in Fig. 1(a) rather than the asymmetric one 
caused by depopulation of -1| sub-band only. Secondly, hole-doping enlarges the distance 
between the Fermi points in the -1J, sub-band as shown in Fig. 1(b), which requires a large 
field to detouch from the Fermi level. On the contrary, our experimental results show smaller 
transition fields for a' in the neutron-irradiated graphite that has a larger hole concentration. 
Thirdly the nearly simultaneous depopulation of two sub-bands seems to be natural in the 
point of view of the charge neutrality condition. Therefore, we interpret the behavior of Hall 
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Fig. 6. (Color online) The proposed sub-band structure realized above 53 T and possible three nesting vectors, 
labeled I, II, and III. 


resistivity and magnetization as features of the realization of the quantum limit state in which 
only two sub-bands (electron-like OJ, and hole-like -It sub-bands) cross the Fermi level. 

The observed enhancement of p zz above 53 T suggests the emergence of a novel density 
wave state between 53 T and 75 T in the quantum limit state. 14 In this state, the possible 
nesting vectors are shown in Fig. 6 by arrows labeled I, II, and III. Here, we suppose simul¬ 
taneous nesting between the residual two Fermi points connected by the same wave vector 
for the each state to reproduce the observed insulating behavior, although the nesting vectors 
are not described in Fig. 6 for clarity. The vector I represents the charge density wave phase, 
whereas the inter-sub-band nesting II and III denote the excitonic phase. 4,5,25,26 According to 
the original argument of the excitonic insulator, 26 the present gapped state caused by the insta¬ 
bility of the Fermi surface can be expressed by the BCS-like formulation, and called excitonic 
BCS-like state. As stated in the Ref. 6, the charge density wave phases can be less stable be¬ 
cause of the direct Coulomb interaction. Therefore, gap will be formed at the Fermi energy 
by inter-sub-band nesting with vector II or III, and simultaneously between the residual two 
Fermi points. Taking into account the charge neutrality condition (k 6 = kp 1 ), the length of the 
vector III (2kjl + 8k) is locked to k/cq. We suppose such a commensurate modulation can be 
favorable than the incommensurate one characterized by the vector II. In this commensurate 
state, the carriers in a pair have finite average momentum, i. e. (kj? +8k)-k® = 6k, analogous to 
the FFLO phase in superconductors. 27,28 With increasing magnetic field, the value of 2+ 6k 
remains n/co, whereas the 6k becomes larger. The reported additional anomaly at 75 T might 
be a transition caused in the process of increasing the 6k. The true nature of this transition 
will be clarified through additional experiments in the future. 
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4. Conclusion 

The transport and magnetic properties of graphite were investigated in pulsed magnetic 
fields up to 75 T. The Hall resistivity becomes nearly zero as the applied magnetic field 
approaches 53 T. In addition, an anomaly is observed in differential magnetization at about 
53 T. This transition field is reduced by doping holes. These results support the realization 
of the quantum limit state, where only two sub-bands (spin polarized one electron-like and 
one hole-like sub-band) cross the Fermi level, in the field region over 53 T. Formation of the 
density wave state in this situation can be regarded as the excitonic BCS-like state. 
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